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Abstract. We describe an algorithm for obtaining generators of the unit 
group of the integral group ring ZG of a finite abelian group G. We used our 
implementation in MAGMA of this algorithm to compute the unit groups of ZG 
for G of order up to 110. In particular for those cases we obtained the index of 
the group of Hoechsmann units in the full unit group. At the end of the paper 
we describe an algorithm for the more general problem of finding generators 
of an arithmetic group corresponding to a diagonalizable algebraic group. 



1. Introduction 

Let G be a finite abelian group. For a ring R we let RG be the group ring over 
i?, consisting of sums X^^gg ^g9^ with Og G R. We take i? = Z and consider the 
unit group 

(ZG)* = {it e ZG I there is a w G ZG with vu = 1}. 

Higman ([II]) showed that (ZG)* — ±G x F, where F is a free abelian group. 
Moreover, Ayoub and Ayoub ( 1 ) have established that the rank of F is \{\G\ + 
1+^2 — 21), where ^2 is the number of elements of G of order 2, and I is the number 
of cyclic subgroups of G. 

Hoechsmann ([15]) described a construction of a set of generators of a finite- 
index subgroup of (ZG)*, called the group of constructable units. Regarding this 
construction he wrote "Does this method ever yield all units if n — \G\ is not a 
prime power? The answer seems to be affirmative for all n < 74." In [TS] this 
question is not dealt with any further. When n = 74 it is known that the group of 
constructable units is of index 3 in the full unit group (see [16]). So the question 
remains whether the constructable units generate the full unit group if |G| < 74. 

In this paper we develop algorithms for computing generators of the unit group 
(ZG)*. Using our implementation of these algorithms in the computer algebra 
system MAGMA ([Ij) we have computed the unit groups for all abelian groups of 
order < 110. We found 12 groups G of order less than 74 whose unit group is not 
generated by the Hoechsmann units (namely, the groups of order 40, 48, 60, 63 and 
65). 

In the next section we start by collecting some well-known facts and immediate 
observations concerning lattices, groups, and associative algebras. Also, in the 
second half of the section (Section [^3]) we describe our approach to computing the 
unit group of the maximal order in a cyclotomic field. This is achieved by combining 
a construction by Greither (\12\) of a finite- index subgroup of the unit group, along 
with a Magma program by Fieker for "saturating" a subgroup at a given prime p. 
The latter algorithm and its implementation will be described elsewhere. 

Section [3] contains the main algorithm of this paper, namely an algorithm for 
computing the unit group of an order O in a toral algebra A. The main idea is to 
split A in its simple ideals CiA, where the et are orthogonal primitive idempotents. 
The CiA are number fields with orders ejC So in order to compute their unit groups 
we can use the effective version of the Dirichlet unit theorem (cf. [S], [IH])- The 
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basic step of the algorithm is, given two orthogonal idempotents ei, 62, to obtain 
the unit group of (ei + €2)0 given the unit groups of eiO, i — 1,2. 

In Section |4] we describe our method for obtaining generators of unit groups of 
integral abelian group rings. Its main ingredients are the construction of the unit 
groups of cyclotomic fields, and the algorithm of Section |31 We comment on the 
running times of the implementation of the algorithm in Magma, and we give a 
table containing all abelian groups of orders up to 110, where the constructable (or 
Hoechsmann-) units do not generate the full unit group. For all these groups we 
give the index of the group of constructable units in the full unit group. 

Finally in the last section we indicate an algorithm to obtain generators of the 
arithmetic group corresponding to a connected diagonalizable algebraic group de- 
fined over Q. Again the main ingredient is the algorithm of Section [31 

In our main algorithms and implementations we make essential use of Fieker's 
implementation in Magma of an algorithm by Ge ( 10 ) to obtain a basis of the 
lattice 

{(ai,...,a„)eZ"|<^---<" = l} 
of multiplicative relations of given elements wi, . . . , m„ in a number field. 

Acknowledgement: we are very grateful to Glaus Fieker for giving us early 
access to his Magma programs MultiplicativeGroup and Saturation, without which 
this project could not have been carried out. 

2. Preliminaries 

2.1. Lattices. In this paper we use the term "lattice" for a finitely generated sub- 
group of Z™. A lattice A c has a basis, that is a subset ui, . . . ,Ur such that 
every u e A can uniquely be written as it = X]i=i Q^i^i: with £ Z. The lattice 
A C Z™ is called pure if Z"/A is torsion-free. (See [6], §111.16.) 

Let A C Z™ be a lattice with basis ui, . . . ,Ur- We form the r x m-matrix B 
with rows consisting of the coefficients of the Ui with respect to the standard basis 
of Z™. By computing the the Smith normal form of B we can effectively compute 
the homomorphism t/- : Z™ ^ Z^/A (cf. [H], §8.3). Let T denote the torsion 
submodule of Z^/A. Then ip~^{T) is the smallest pure lattice containing A. So in 
particular, the Smith normal form algorithm gives a method to compute a basis of 
the lattice V f) Z™, where F is a subspace of Q™. For example, we can compute 
the intersection of lattices this way. 

As observed in [5], §111.16, a lattice A is pure if and only if it is a direct summand 
of Z™. So in that case, by computing a Smith normal form, we can compute a basis 
of Z™ such that the first r basis elements form a basis of A. 

Remark 2.1. Gomputing a basis of the lattice V D Z™, where V is a subspace of 
Q™, is called saturation. There are methods known for this, based on computing 
a Hermite normal form, that are more efficient than the approach outlined above. 
Going into these matters would lead us too far from the subject of this paper. 

2.2. Toral algebras. We say that an associative algebra A over Q is toral if it is 
semisimple, abelian and has an identity element, which we will denote by e. For 
example the group algebra QG of a finite abelian group G is toral. 

By the Wedderburn structure theorem (cf. (TS], §3.5) a toral algebra A is a 
direct sum A = Ai ® ■ ■ ■ (B Ag, where the Ai are ideals that are isomorphic (as 
associative algebras) to field extensions of Q. A nonzero element Cq G A is said 
to be an idempotent if e§ = eg. Two idempotents 61,62 are called orthogonal if 
6162 = 0. Furthermore, an idempotent is called primitive if it is not the sum of 
orthogonal idempotents. Now the decomposition of A into a direct sum of simple 
ideals corresponds to a decomposition of the identity element e G A as a. sum of 
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primitive orthogonal idempotents, e — ei + ■ ■ ■ + Bs- Here is the identity element 
of Ai, and vice versa, Ai = eiA. We remark that there are algorithms to compute 
the Ci, given a basis of A (cf. [8 , [9 ). 

When A — QG, with G a finite abelian group, there is a very efficient way to 
compute the primitive idempotents. Let x : G — > C* be an irreducible character of 
G. Then 

= 1777 XI 

' ' gee 

is an idempotent in CG. Moreover, the e^^, as x runs over all irreducible charac- 
ters, are primitive orthogonal idempotents with sum e (cf. [B], Theorem 33.8). In 
particular, they form a basis of CG. Let m denote the exponent of G, then all 
irreducible characters x have values in the cyclotomic field Q(Cm)- So the Galois 
group Gal((Q(<;^m)/Q) = (Z/mZ)* acts on the irreducible characters. Now we sum 
the e^, for x in orbit of Gal(Q(Cm)/Q), and obtain the primitive orthogonal 
idempotents of QG. 

A subset O of a toral algebra A, containing the identity of A, is said to be an 
order (or, more precisely, a Z-order) if there is a basis ai, . . . , am of A such that 
O ~ Zai + • • • + Zcm and aia^ G O for 1 < i, j < m. For example, ZG is an order 
in QG. The unit group of O is 

O* = {a e O I there is 5 € O with ab = e}. 

We consider the problem of obtaining a basis of the lattice of multiplicative relations 

L = {{ai, . . . , a^) e Z'- I ar • • • = e}, 

where ai, . . . ,0^ are given elements of O* . Note that CiA are number fields. So 
using Ge's algorithm (fTU]) we can compute bases of the lattices 

Lj = {(ai, . . . , a^) e Z'' I (ej-ai)"! ■ • ■ (ej-a^)"'- = e^}. 

Moreover, L = HjLj so we can compute a basis of L (see Section [^?T|) . 

2.3. Standard generating sets. Let J7 be a finitely-generated abelian group. We 
say that a set of generators 171 , . . . , of J7 is standard if 

(1) for 1 < i < s the order of gi is di, 

(2) for s + 1 < i < r the order of is infinite, 

(3) di\d,+i for i < s, 

and there are no other relations. So a standard set of generators immediately gives 
an isomorphism of U to Z/diZ © • • • ® Z/d^Z ® Z'^". 

Giving finitely-generated abelian groups by standard generating sets yields straight- 
forward algorithms for several computational tasks concerning these groups, such as 
computing the index of a subgroup, and computing the kernel of a homomorphism. 

If an abelian group is given by a non-standard set of generators, then we can 
compute a standard one by computing the lattice of all relations of the generators, 
followed by a Smith normal form computation (cf. [21^, §8.3). So, using the algo- 
rithm indicated in the previous section, we can compute a standard set of generators 
for a finitely-generated subgroup of O*, where O is an order in a toral algebra. 

Remark 2.2. For many computational problems regarding finitely-generated abelian 
groups it suffices to compute a Hermite normal form of the relation lattice. How- 
ever, in our applications the main computational problem is to obtain the relation 
lattice (see Tabled]), the subsequent computation of the Smith normal form does 
not bear heavily on the running time. Therefore, for our purposes, a Smith normal 
form is the most convenient. 



4 



PAOLO FACCIN, WILLEM A. DE GRAAF, AND WILHELM PLESKEN 



2.4. Units of cyclotomic fields. Let n be a positive integer. We consider the 
cyclotomic field Q(Cn), where Cn is a primitive 71-th root of unity. The ring of 
integers of this field is Z[C„]. By Dirichlet's unit theorem the unit group 
is equal to T x F, where is a free abelian group of rank |(/?(n) — 1, and T is 
the group of roots of unity in Q(Cn)- The problem considered in this section is 
to obtain generators of the unit group, There are algorithms for this that 

work for any number field (cf. [S], |19)): but their complexity is such that it is only 
practical to use them for n up to about 20 (depending on the hardware one uses, 
of course). For this reason we sketch a different approach, using several results 
from the literature. The situation is straightforward when n is a prime power, see 
Section 12.4.11 Then in Section 12.4.21 we describe what can be done when n is not 
a prime power. Using these methods we obtained a list of the generators of the 
unit groups Z[(^„]*, for n < 130. However, for several n the correctness of this list 
depends on the Generalised Riemann Hypothesis, that is for n prime between 67 
and 127, and for n = 115, 119, 121, 123, 125, 129 (so 19 cases in total). 

Throughout we set Q(C„)+ = ]RnQ(C„); then Q(C„)+ = Q(C« + C')- By h+ we 
denote the class number of Q(Cn)^- 

2.4.1. When n is a prime power. Suppose that n = is a prime power. For 
1 < a < f with gcd(a,p) — 1 set 

1 — g 1 /^a 

C — /- 2 ^ 

Then hes in the unit group of Q(Cn)^- Let t/„ be the group generated by — 1, 
(^n and all £,a- Then for the index we have : C/„] = /i+ (this is obtained 

by combining Corollary 4.13, Lemma 8.1 and Theorem 8.2 in [^). It is known 
that — 1 a f{n) < 66, and assuming the Generalised Riemann Hypothesis, we 
have = 1 when if{n) < 162 (see the Appendix in 22J). So for those n we have 
generators of the unit group. 

2.4.2. When n is not a prime power. Here the situation is more difficult. First of 
all we assume that n ^ 2 mod 4, as for n = 2 mod 4 we have that Q(Cn) and Q(Cn) 

are isomorphic. By £'+ denote the unit group of Q(Cn)''"- We use a finite-index 
subgroup oiE'^ defined by Greither ([H]). Here we briefly describe his construction. 

Let Q = Gal(Q(Cn)/Q), and write the elements of Q as Ua, where gcd(a,n) = 1 
and aa(Cn) = Cn- For a^Y^a ^aCTa e 'ZG and x e Q(Cn) define 

x'' = l[aa{xr^. 

a 

Let n = Yli=i pT t>e the factorisation of n in distinct prime powers. Set S = 
{1, . . . , s} and Ps = {/ C 5 I / 7^ S}. For I € Ps we set m = UieiPT- 

We consider arbitrary maps /? : — )• ZCJ, which we extend to maps (denoted by 
the same symbol) l3:Ps^Zghy /3(0) = 1, /3({i}) = l3{i), and /3(/U J) = P[I)P{J) 
if /nJ = 0. Now let z e Q(C„). For / e Pg set z/ = 1-z"^ and z(/3) = [l/eS^f 
Set t — — X]/es "^/^(^) ^ Then for a with \ < a < ^ and gcd(a,n) = 1 we 
consider 

ia{P) = Cn ^^^^ : ^hcrC da = . 

1 

(Note that for n odd, Cn lies in Z[C„], whereas for n even we have a odd and hence 
is an integer.) 

Following Greither we describe a good choice for /3. First a small piece of nota- 
tion: if g is an element of order m of a group, then we set Ng = 1 + g + ■ ■ ■ + g™'^^, 
which lies in the corresponding integral group ring. Now consider an i d S. Let Gi 
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denote the Galois group of Q(Cn/p''0^ over Q. This group contains the Frobenius 
automorphism Fi (by definition: Fi{(^^p^i) — C'^'/p'i)- This yields the element Np. 
in "EQi. Now we define /3(z) to be a lift of Np. to ZQ. 

Let Cp be the subgroup of generated by —1 and the Ca(/3)- Greither proved 
that Cp does not depend on the choice for the lifts of the Np., and that it is of 
index h'^ip in , where 

^p = f[er'ff'^'\ 

1=1 

here e^, fi and gi are respectively the ramification, inertial and decomposition 
degree of in Q(Cn)"'" (so that Cifigi = ^(p{n)). 

Now let Un denote the group generated by (n and Cp . Then using [22] , Corollary 
4.13, we get that [Z[C„]* : Un] = 2h+ip. 

We used a program that Glaus Fieker has written in Magma, V2.17-2 (see [2] for 
some of the background). This program, given a finite-index subgroup V of Z[Cn]*, 
and a prime p, computes a subgroup V of Z[C„]*, containing and such that the 
index [Z[(^„]* : V] is not divisible by p. We used this starting with the group Un- For 
n < 130 and (p{n) < 72 we have hn = 1, and assuming the Generalised Riemann 
Hypothesis, we have /i+ = 1 for all n < 130 ([12, Appendix). Therefore we can 
compute [Z[<^„]* : [/„], and get all primes dividing it. So in the end we arrived at 
the full unit group Z[C„]* for all n < 130. 

3. Unit groups of orders in toral matrix algebras 

Let A be a toral algebra with identity e, and O C A an order. In this section 
we consider the problem of computing a set of generators of O*. First we consider 
some special cases. 

3.1. A simple toral algebra. Let A be a simple toral algebra. In other words, it 
is isomorphic to a finite extension of Q. Now there are algorithms for computing 
generators of an order in a number field (the effective version of the Dirichlet unit 
theorem, see [S], [T^). We use these to compute generators of O* as well. 

3.2. Two idempotents. Let A be a toral algebra with identity e, and ei,e2 S A 
orthogonal (but not necessarily primitive) idempotents with ei + 62 — e. Set 
Ai = BiA, then A = Ai (B A2. Let O be an order in A, then Oi = CiO is an order 
in Ai. Here we suppose that we have generators of O*, i = 1, 2, and the problem is 
to find generators of O*. 

Set J — {eiO n O) + (62© n O). Since this is an ideal in O we can form the 
quotient R — O/J. Consider the maps ipi : e^O — R defined by (pi{eia) = a + J, 
for a E O. (Note that this is well-defined: if eia — eib then a — b = 62(0 — b) 
so it lies in J.) These are surjective ring homomorphisms with respective kernels 

aO n O. 

Lemma 3.1. O = {ai + a2 \ ai G CiO and (pi{ai) — ip2{a2)}- 

Proof. Let a G O and set Oi — Cia then a — ai + a2 and G e^C Moreover 
tpi{ai) = a + J = ip2{a2). Conversely, let a, 6 G O be such that (^i(ai) ~ f2{a2), 
where oi = eio and 02 = 626. Then a — b € J, whence a — b + ui + U2 where 
Ui G CiOnO. Therefore eio — eib+ui so that ai+a2 = eia + e2b — eib+e2b+ui = 
b + ui which lies in O. □ 

Corollary 3.2. O* = {ai + 02 | G {ciO)* and fi{ai) — <^2(a2)}- 

So in order to compute generators of O* we perform the following steps: 
(1) Compute bases of Oi = CiO. 
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(2) Compute bases of Oi n O, of J = (d n O) + {O2 n O), and set B = O/J. 

(3) Compute generators of the groups Hi ~ ipi{0*) C R* and H = Hi Ci H2. 

(4) Compute generators of the groups Mi = ip~^{H) C O*. 

(5) Compute generators of the group O* = {ai + 02 | G Mi and ipi{ai) = 
■^2(02)}- 

3.2.1. Implementation. We comment on the implementation of the steps of the al- 
gorithm. Step (1) is done by a Hermite normal form computation. The intersections 
in Step (2) are computed using the techniques indicated in Section \TJ\ Note that 
a basis of J is obtained by concatenating the bases of Oi fl O. The ring R = O/J 
can be constructed by a Smith normal form computation. 

For Step (3) we assume that is isomorphic to a number field (in other words, 
that 62 is a primitive idempotent). When using the algorithm, this can always be 
arranged (see Section . Then 62© an order in it. We set / = 62© fl O and use 
the isomorphism R = (62 0)//. Subsequently we use algorithms described in |13j . 
[T7] to compute a standard generating set of {e20/I)* . So computing Hi, H2 as 
subgroups of (620//)* we can perform the operations of Step (3). 

In Step (4) we view ipi as a homomorphism O* — >■ Hi. We use this to compute 
generators of the kernel of (.pi as well as pre-images of the generators of H . Together 
these generate the group Mi. As remarked in Section [531 '^e can compute standard 
generating sets of the groups O*. Using these, it is straightforward to obtain a 
standard generating set for the subgroups Mi. Then we restrict (pi to obtain a 
homomorphism Lpi : Mi H . 

Now we come to Step (5). Let hi,...,hr, ai, . . . ,as, bi,...,bt be standard 
generating sets of H, Mi and M2 respectively. Set 

A = { (ai , . . . , a„ /3i , . . . , A) e | ^1 (a^ • • • = ^2 (6^ • • • ) }• 

Then 

O* = {a^ • • • + 5f ^ • • • 6fM ■ • • , • • ■ , A) e A}. 

Moreover, A is a lattice, hence has a finite basis. Furthermore, the elements of 
O* corresponding to the elements of a basis of A generate O* . So the problem of 
finding a generating set of O* is reduced to finding a basis of A. 
Define iiij , i>ij 6 Z by 

r 

^i(ao=n'^r 

r 

i=i 

A small calculation shows that (ai, . . . , as, . . . , /3t) £ A if and only if 

s t 

(1) ^ ^J.ija^ - ^ VkjlSk = mod ord(/ij) 

i=l fc=l 

for 1 < j < r (and where ord(hj) denotes the order of hj). Let S be the integral 
matrix with columns 

{fiij, Hs3,-i^ij, ■ • • , -iytj,ord{hj)) 

for 1 < j < r. We compute a basis the integral kernel of S, which is the set of all 
V e Z*+*+^ such that vS = 0. For each v in this basis we take the vector consisting 
of the first s + t coordinates. This way we obtain a basis of A. 
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3.3. The general case. Now let yl be a toral algebra, ei,...,em its primitive 
orthogonal idempotents with sum e, and Ai — eiA the corresponding simple ideals. 
Let O be an order in A] then CiO is an order in A^, and as indicated in Section 
13.11 we can compute generators of the unit groups (ciO)*. Then for j — 2,3, .. . 
we set Cj = ei + ■ ■ ■ + ej and we apply the algorithm of Section 13.21 to the algebra 
tjA with its order ejO, and two idempotents ^j-i and Cj, yielding the unit group 
(tjO)* . When the algorithm terminates we have the unit group O* . 

4. Units of integral abelian group rings 

Let G be an abelian group. As seen in Section 12. 2[ it is straightforward to 
compute primitive orthogonal idempotents ei,...,er G QG such that ei(QG) is 
isomorphic to a field extension of Q. Let m denote the exponent of G, then (cf. [1]) 

td 

d\m i—1 

where Q(Cd) is the cyclotomic field of order d, and td is the number of cyclic 
subgroups of G of order d. In particular, ei(QG) = Q(Cdi)- We consider the order 
ZG in QG. We have that ei(ZG) is isomorphic to ^[Cd.], and by the results of 
Section we have generators of Z[^rf.]* for di < 130. So for small groups G we 
can apply the algorithm of Section [3] to obtain generators of the unit group (ZG)*. 
Using our implementation of the algorithms in Magma, we have carried this out 
for all abelian groups of orders up to 110. In Table 2] we collect some timings and 
other data related to the algorithm. 



G 


(^(exp(G)) 


1 digits] 




^tot 


G70 


24 


7.4 


314 


328 


Gso 


32 


63.4 


1131 


1183 


G90 


24 


159.5 


1043 


1078 


G91 


72 


3.7 


2352 


2446 


Gge 


32 


31.2 


2322 


2373 


G2 X G48 


16 


181.3 


1575 


1781 


G2 X G2 X G24 


8 


54.6 


1031 


1267 


G2 X G4 X G12 


4 


22.5 


537 


725 


Gioo 


40 


217.6 


3822 


3942 


G2 X G50 


20 


66352.6 


414569 


426654 


G5 X G20 


8 


379.7 


1227 


1425 


Gio X Gio 


4 


275.2 


1131 


1332 



Table 1. Runtimes (in seconds) for the algorithm to compute 
generators of (ZG)*. The first column lists the isomorphism type 
of G, and the second the value of the Euler (/j-function on the 
exponent of G. The third column has the average number of digits 
of the coefficients of the units output by the algorithm (with respect 
to the standard basis of ZG). The fourth column displays the time 
spent to compute multiplicative relations in cyclotomic fields. The 
last column has the total time spent by the algorithm. 

From Table [3] we see that the running time is dominated by the time needed to 
compute multiplicative relations in cyclotomic fields (Ge's algorithm). This algo- 
rithm needs to work harder if the degrees of the fields that occur are higher. Indeed, 
the running times generally increase when (p{ex-p{G)) increases (note that this is 
the highest degree of a cyclotomic field occurring in the decomposition of QG). 
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However, also the size of the elements of which we need to compute multiplicative 
relations plays a role. For some groups the average number of digits of a unit, as 
output by the algorithm, is very high. This is seen most dramatically for C2 x C50. 
Note also that the size (i.e., the average number of digits of their coefficients) of the 
units output by the algorithm is far from being optimal; indeed, for G = C2 x C50, 
the unit group (ZG)* is also generated by the Hoechsmann units (see below). 

Various constructions of finite index subgroups of (ZG)* have appeared in the 
literature (see [3D] for an overview). Among these a construction by Hoechsmann 
([15)1 seems to yield subgroups of particularly small index. As an application of our 
algorithm, we have compared the full unit groups with the groups of Hoechsmann 
units, for G of size up to 110. 

We briefly describe Hoechsmann's construction. First, let G be a cyclic group of 
order n, generated by the element x. For i>{) and y G C we set 

Let i,j be integers with < i, j < n and gcd(i,rt) — gcd(j, n) — 1. Let k,l be 
positive integers with li = 1 + kn. Set 

Then Uij{x) is a unit in (ZG)*. Let 6(G) denote the set of all units constructed 
in this way. We let H be the group of units in (ZG)* generated by all 0(G), where 
G ranges over the cyclic subgroups of G of order > 2, along with ±G. It is called 
the group of constructable units of ZG. By 15j, Theorem 2.5, H is a subgroup of 
finite index of (ZG)*. 

In general the set consisting of the Uij{x) is a heavily redundant generating set. 
We note, however, that [H] Theorem 1.1, gives a non- redundant generating set of 
H, if the group Hm = (Z/mZ)*/{±l} is cyclic, where m is the exponent of G. We 
remark that for small m the group Hm is cyclic with only few exceptions. The 
values of m up to 120, for which Hm is not cyclic, are 24, 40, 48, 56, 60, 63, 65, 72, 
80, 84, 85, 88, 91, 96, 104, 105, 112, 117, 120. 

We let Hind(G) be the index of H in (ZG)*, and we call this number the Hoechs- 
mann index. We have computed the Hoechsmann indices for all abelian groups of 
orders up to 110. For most groups the index is 1. The groups for which it is not 1 
are listed in Table |4l along with the corresponding Hoechsmann indices. 

From Table |4] we see that on many occasions when Hind(G) ^ 1 we have that 
|G| = m, with Hfn not cyclic. Also, for the groups considered, if for one group G 
of order m we have Hind(G) 7^ 1, then the same holds for all groups of that order 
(except when (ZG)* has rank 0). 

Remark 4.1. For the groups Gp, with p a prime between 67 and 120 the correctness 
of our computation depends on the Generalised Riemann Hypothesis (see Section 

[23). 

5. Computing generators of arithmetic groups of diagonalisable 

algebraic groups 

In this section we use the term D- group as short for for diagonalisable algebraic 
subgroup of GL(n, C), defined over Q. We consider the problem to compute gener- 
ators of the arithmetic group G(Z) = G fl GL(7t., Z) for a given connected D-group 
G. We assume that the group G is given by its Lie algebra g C gl{n, C) which in 
turn is given by a basis of matrices with coefficients in Q. In the sequel such a basis 
will be called a Q-basis. 

We have two subsections: the first is devoted to characters, and in the second 
we outline our algorithm. 
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G 


ffind(G) 


G 


Hind(G') 


G40 






2 


Gsi 






2 


C2 X 


(^20 




2 


G2 X 


C42 




2 


C2 X 


C2 X 


Gia 


2 








2 


C48 






2 


C90 






3 


C2 X 


C24 




2 


C3 X 


C30 




3 


C2 X 


G2 X 


G12 


2 


Ggi 






3 


C4 X 


C12 




4 


C96 






8 


Geo 






2 


C2 X 


C48 




16 


G2 X 


C30 




2 


Ca X 


C2 X 


^24 


32 








3 


G2 X 


G2 X 


G2 X C12 


16 


C3 X 


G21 




3 


C2 X 


C4 X 


^12 


64 








2 


Gi X 


C24 




64 


C74 






3 


C98 






7 


Cso 






4 


C7 X 


Ci4 




343 


C2 X 


GiQ 




8 


^104 






2 


C2 X 


G2 X 


G213 


16 


G2 X 


^52 




2 


C2 X 


C2 X 


G2 X Cio 


32 


C2 X 


C2 X 


G2Q 


2 


C4 X 


G20 




8 


Ci05 






4 



Table 2. Hoechsmann indices for abelian groups of orders up to 
110, for which this index is not 1. 



5.1. Computing the characters. We recah that a character of a group ff is a 
homomorphism x '■ H ^ F* , where F* denotes the multipUcative group of nonzero 
elements of a field F. In this section we let G C i? C GL„i(C) be two connected 
D-groups defined over Q. A well-known theorem (cf. |3j, Proposition 8.2) states 
that there are characters Xi; ■ • ■ jX* ■ ^ ~^ *C* such that G is the intersection of 
their kernels, i.e., 

G ^{heH \ Xi{h) = 1 for 1 < i < s}. 

The problem considered in this section is to find a number field if, and characters 
Xi, . . . ,Xt • ~^ K* such that G is the intersection of their kernels. We require 
that K and the Xi be explicitly given, i.e., that we know the minimum polynomial 
of a primitive element of and that for a given h € H we can effectively compute 
Xi{h)- We assume that we are given the Lie algebras i),g oi H and G respectively, 
in terms of Q-bases. 

Proposition 5.1. Let X E GL„j(C) be such that H = XHX^^ consists of diagonal 
matrices. Set t) — Xt)X^^, g ~ XgX^^. Then consist of diagonal matrices as 
well. Set 

Hh) = {(ei, ■ • ■ ,em) e Z™ I eiai H h e„ia,„ = for all diag(Q;i, . . . ,a„i) £ t)}, 

and similarly define A(§). Then A(f)) C A(g) C are pure lattices (cf. Section 
\2.1\) . Let e^ = (e;^!, . . . , 6^^^) e Z™, for 1 < i < s, be a basis of A{q) such that 
Cj+i , . . . , is a basis 0/ A(()). For 1 < i < t define Xi ■ H ^ C* by XiW = 
c^i''^ •••am™, where XhX^^ — diag(Q!i, . . . , a„i)- Then G, as subgroup of H, is 
the intersection of the kernels of the Xi ■ 

Proof. Note that [} is the Lie algebra of H . Therefore, it consists of diagonal 
matrices. Since g C f) we get the same property for g. The lattices A(t)), A(g) are 
obviously pure. So as noted in Section [^Tl the required basis of A(g) exists. 

A character x of -ff is given by a sequence (ei,...,e„i) of integers such that 
x(diag(ai, . . . , a„i)) = cii ■ ■ ■ cttn ■ This follows from the same statement for the 
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full diagonal group Dm, along with the fact that a character of H extends to 
a character of D„i ([5], Proposition 8.2). Furthermore, the differential of x is 
given by dx(diag(ai, . . . ,am)) = eiai + • • ■ + emCtm- So since f) consists of all 
diag(Q;i, . . . , am) with eiai + • • • + CmCtm = for all (ei, . . . , e™) € A(f)) we get 
that 

H = {diag(ai, . . . ,a™) | a"j^ ■ • • =1 for aU (ei,...,em) € A(^)}. 

Indeed, the latter is a connected algebraic group with the same Lie algebra as H. 
An analogous statement holds for G — XGX~^ and g. This implies the statement 
of the proposition. □ 

So in order to solve the problem of this section we need to say how the various 
objects in Proposition 15. II can be computed. 

Let B C AIm{Q) be the associative algebra with one generated by the basis 
elements of () . Then, using the terminology of Section 12.21 B is toral (indeed: its 
elements can simultaneously be diagonalized) . Now using algorithms given in for 
example [7], we can compute a splitting element bo G B. This means that 
bo generates B (as associative algebra). By repeatedly factoring polynomials over 
number fields, we construct the splitting field K of the minimal polynomial of 
bo over Q. After factoring this polynomial over K we find the eigenvalues and 
eigenvectors of bo- From this we get a matrix X € GLm(K) such that XboX~^ is 
diagonal. We note that this also implies that Xi)X^^ consists of diagonal matrices. 
As this is the Lie algebra oi H — XHX^^, we get the same for H, as the latter 
group is connected. 

Since we have given a basis of f) we can compute a basis of X\jX^^. Cor- 
responding to a basis element diag(Q;i, . . . ,am) we construct the linear equation 
aiXi + ■ ■ ■ amXm = in the unknowns Xi. We can express the as vectors in Q"^, 
where d is the degree of K. So by solving these equations we get a basis of the 
space 

V{[)) = {(ai, . . .,am) e I aiai H \- amam = for aU diag(ai, . . . ,a„i) G [)}. 

Now A{[)) = V {i)) n Z™, so as seen in Section [23 we can compute a basis of A{i)). 
Furthermore, using the methods outlined in the same section, we get a basis of A(g) 
that contains a basis of A(f)). Then the characters Xi ^re easily constructed as in 
Proposition l5.ll 

5.2. Putting the pieces together. Now let G C GLm(C) be a connected D- 
group, defined over Q. We assume that we have given its Lie algebra g C 0(,„(C) 
by a Q-basis. We describe an algorithm for obtaining a finite set of generators of 
the group G(Z) = G n GL„(Z). 

Lemma 5.2. Let A C Afm(C) be the associative algebra with one generated by the 
basis elements of g. Let A* denote the set of invertible elements of A. Then A* is 
a D-group defined over Q with Lie algebra A. Furthermore, G G A* . 

Proof. Note that A has a basis consisting of matrices with coefficients in Q. Let 
A{Q) denote the Q-algebra spanned by such a basis. Consider the linear equations 
that define A{Q) as a subspace of Mm(Q). Then those equations define A* as 
a subgroup of GLm(C). Since the basis elements of g can be simultaneously be 
diagonalised, the same holds for the elements of A*. It follows that A* is a D- 
group, defined over Q. 

Since A* , as subgroup of GLm(C), is given by linear equations, its Lie algebra, 
as subalgebra of g(^(C), is given by the same equations. Hence the Lie algebra of 
A* is A, where the Lie product is given by the commutator. So since g C A, we 
also have G C A*, as G is connected. □ 
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Now in order to conipute generators of G'(Z) we have the following algorithm. 
As input we take a Q-basis of g. 

(1) Compute a Q-basis of A (notation of Lemma r5.2[) . 

(2) Let O = An Mm(Z) and use the algorithm of Section [3] to compute gener- 
ators ai , . . . , Or of O* . 

(3) Use the algorithm of Section 15.11 to construct a number field K and char- 
acters xi, . . . , Xs : A* K* such that G is given, as subgroup of A* , by 
the intersection of their kernels. 

(4) Use Ge's algorithm ([IH]) to compute bases of the lattices 

A,-{(ei,...,e,) eZ'- |x^(air---X.Kr- = 1}. 

(5) Let ej, — {ek,i, ■ ■ ■ ,ek.r), for 1 < fc < i, be a basis of the intersection of 
all Ai (1 < i < s). Set Uk = al'''^ ■ ■ ■ a'r''''' . Then ui,...,ut is a set of 
generators of G(Z). 

Proposition 5.3. The previous algorithm terminates correctly. 

Proof. Note that by Lemma [5.21 A is the Lie algebra of A*, so that in Step (3) 
we have the correct input for the algorithm of Section 15.11 A basis of O can be 
computed using the methods of Section 12.11 So all steps can effectively be carried 
out. 

Observe that A*{Z) ^ A* n GL,„(Z) is equal to O* . Now g e G(Z) if and only if 
g G A*{'L) and xAq) — 1 for 1 < i < s. But this is equivalent to g = ■ ■ ■ a'^^ for 
certain £ Z, and Xi(g) = 1 for 1 < i < s. But the latter condition is equivalent 
to (ei, . . . , Cr) e Ai for 1 < i < s. Hence g G G(Z) if and only if g is a product of 
the Uk. □ 
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